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Abstract 

The quantum dimer model, relevant for short-range resonant va¬ 
lence bond physics, is rigorously shown to have long range order in a 
crystalline phase in the attractive case at low temperature and not too 
large flipping term. This term flips horizontal dimer pairs to vertical 
pairs (and vice versa) and is responsible for the word ‘quantum’ in 
the title. In addition to the dimers, monomers are also allowed. The 
mathematical method used is ‘reflection positivity’. The model and 
proof can easily be generalized to dimers or plaquettes in 3-dimensions. 


1 Introduction 

Starting with the work of Rokhsar and Kivelson d! quantum dimer models 
became popular in connection with the study of resonating valence bond 
(RYB) states, as well as in the study of cold bosons, frustrated magnetism, 
Josephson junction arrays, and other physical models, see [13] for a review. 
In these models dimers on a square lattice not only have a pair interaction 
that favors or disfavors parallel dimers next to each other, but they also have 
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a flipping term that exchanges a horizontal pair with a vertical pair. It is 
this ‘non-diagonal’ flipping term that gives rise to the appellation ‘quantum’. 

Many kinds of physical phases can arise from such models. As the in¬ 
teraction is increased from large negative values (strongly attractive case) to 
large positive values (strongly repulsive case), quantum dimers are expected 
to exhibit a remarkable series of transitions: for instance, on the square lat¬ 
tice, the system is expected to pass from a crystalline “columnar” state, to 
a “mixed” state, then a “plaquette”, and finally a “staggered” state, pass¬ 
ing through an anomalous liquid phase at the Rokshar-Kivelson point, see 
m [I6j for details about these phases and transitions. Up to now, none 
of these conjectured phases and transitions has been rigorously established. 
For recent rigorous results about the anomalous liquid phase in the classical 
interacting dimer model, see m 

Here we study the attractive case on the square lattice and rigorously 
establish that there is a columnar state at low temperature provided the 
monomer density and flipping rate are not too large. Our analysis easily 
extends to the three-dimensional cubic lattice, in which case the dimers can 
remain as dimers or they can be replaced by two-dimensional plaquettes. 

Much before Ref. H7j, Heilmann and Praestgaard [10, [TTj, 12], using a 
Peierls contour argument, proved that the same model, but without the 
flipping term, has columnar long-range order (LRO) at low temperature and 
not too large monomer density. 

The introduction of the quantum feature complicates the proof of LRO 
considerably. Borgs et al and Datta et al DEI proved, quite generally, 
that quantum perturbations of classical models do not destroy LRO at low 
temperature. Their methods are both based on a quantum extension of 
Pirogov-Sinai theory [13], 

In our case the method we use - reflection positivity (RP) and chessboard 
estimates - is much simpler. While it was introduced in quantum held theory 
by Osterwalder and Schrader [L|j, its first application in classical statistical 
mechanics goes back to Frohlich, Simon and Spencer [7], and to Dyson, Lieb 
and Simon for the quantum case [3]. See also [6j for an application to the 
anisotropic quantum Heisenberg models, which is relevant for our proof here. 
The use of RP has added advantages apart from brevity. When applicable 
it probably gives better estimates on the range in which LRO holds. In 
(12.21) we give an explicit sufficient condition for LRO. It also establishes 
some correlation positivity results, which can be used for some correlation 
inequalities. 
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2 The model 


We introduce a quantum dimer model on the 2D square lattice (which can be 
easily extended to 3D): consider a square portion of Z 2 of even side L with 
periodic boundary conditions (for convenience later we will choose L to be 
divisible by 4). The Hamiltonian is conventionally expressed in the following 
way: 


H\ = —e (Vphp + h} p vp) — n (VpVp + h^php) + zM = —eT — jiU + zM, 

PC A PcA 

(2- 1 ) 

where the summations run over the plaquettes P of A. The physical objects 
on the lattice are dimers, which lie on the edges, and the symbols hp (resp. 
Vp) denote operators that kill two horizontal (resp. vertical) dimers occupy¬ 
ing P. Correspondingly, hp and v' P create pairs of parallel dimers. Of course 
P can have no dimers, or one dimer, in which case hp and Vp annihilate the 
corresponding state; this means that hp\4>) = Vp\<j>) = 0 if \(j>) is a classical 
dimer configuration such that P has zero or one dimers. We do not assume 
complete dimer covering: if N v and N^ denote the number of dimers in the 
vertical and horizontal directions, M — |A| — 2 N v — 2 Nh denotes the number 
of monomers (empty sites). 

The first term in (12. ip switches the dimers in a doubly occupied pla- 
quette, while the second simply counts the number of such plaquettes. 3D 
extensions of the model, treatable by our methods, are, e.g., the following: 
(i) same Hamiltonian in 3D, where now plaquettes can be in three possible 
directions; (ii) same thing, but with the interaction term replaced by a simi¬ 
lar one, counting the number of unit cubes completely filled by four parallel 
dimers; (iii) quantum Hamiltonian for plaquettes (rather than dimers): par¬ 
allel neighboring plaquettes interact and can flip whenever they occupy the 
same unit cube. See section [3j 

The sign of the constant e is conventional, because it can be changed by 
the unitary transformation i Nv . Therefore, we will assume that e > 0. The 
sign of the second term is very important: we shall be concerned with the 
case that fi> 0. By rescaling, we assume \i — 1. In these units, we assume 
that the monomer chemical potential, z, is larger than —1/2. 

The Hamiltonian preserves the number of dimers (or monomers) and this 
divides the Hilbert space into sectors, each of which has a spectrum. At 
e = 0, there are four ground states, called columnar states, and these are 
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in the sector where the number of dimers is |A|/2. One of them consists of 
parallel columns of horizontal dimers, and the other three are obtained by 
shifting the columns one lattice space to the right, or else rotating the whole 
picture(s) by 90°. 

As anticipated above, we prove that columnar LRO, in the sense defined 
precisely in (I2.10p . is present at small enough e > 0 and large enough /3, 
by using quantum reflection positivity, chessboard estimates, exponential 
localization and an adaptation of the classical Peierls argument. 

More exactly, we prove that a sufficient condition for LRO is that the 
following inequality is satisfied by the three parameters f3, £ and z\ 

' S ^2 / n 2 a n < ^ with a = 18 max -) 24 , e ~^( z+1 / 2 ) j 

n>2 ' 

( 2 . 2 ) 


2.1 Reflection positivity 

We first show that the Hamiltonian (12. ip in a suitable representation is reflec¬ 
tion positive. We associate with each point x e A a Hilbert space "H x = C 5 ; 
there are 5 allowed states at x, one of them is the monomer, and the other 
4 refer to the half-dimers pointing in the directions north, east, south, west. 
The total Hilbert space is the tensor product of these local spaces. To ensure 
that we do not have unpaired half-dimers, we will add to our Hamiltonian a 
term that suppresses unpaired half-dimers. At each point we introduce the 
operators if, with # 6 {0, N, E, S, W}, which are projection operators onto 
the 5 different states, the rotation operator , which flips the state N 

into E and gives zero otherwise, and similarly R E ^ W , etc. 

Consider the Hamiltonian: 


m = a 


iAi-2y(iji^ i + i~q +& ) 


E 1 "! +=E X ° (2.3) 


v + h.c.) (2.4) 


P=(X1,X2,X3,X4) 

V (t N t N t s t s +t E t w t w t E ) 

/ S V X1 JL X2 X3 JL X4 1 X1 J1 X2 X3 JL X4/ 

P=(xi,x 2 ,x 3 ,x 4 ) 


(2.5) 


where the sums in the last two lines runs over the plaquettes of A, each 
plaquette P being thought as the union of the four sites xi,x 2 ,x 3 ,X 4 at the 
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left bottom, right bottom, right top, left top of P, respectively. As we let A —> 
oo, H™ gives infinite energy to unphysical configurations (i.e., configurations 
with unpaired half-dimers) and on the physical space it becomes equivalent 
to (12.ip at /x = 1: note, in fact, that the term in square brackets on the first 
line is zero on physical configurations and positive otherwise; the remaining 
terms (12.41) and (12.51) reduce on the physical space to the first and second 
terms in (12.11) . and similarly reduces to the last term in 02.ip . 

In this representation, there is a natural notion of reflection positivity. 
Take a pair of planes passing through the bonds and cutting A in two equal 
sized halves, A_ and A + . Reflection about these planes consists in changing 
the site labels into their reflected images, and exchanging up with down (N 
with S) or right with left (E with W), depending on whether the planes are 
horizontal or vertical. We call 6 such a reflection; it satisfies the conditions 
listed in [5J. Given any pair of reflection planes, we can write as A + 
9A + Ci9Ci, with A representing the restriction of —H\ to one of the 
two halves, say A + , and CiOCi being the terms connecting the two halves. 
Therefore, H\ is Reflection Positive (RP), which follows from the general 
theory in [[5], Sections 2 and 3.3]. This means that if we take any operator F 
supported on A + then 

{F\9F)) x > 0, (2.6) 

where ((-)) A = Trwith Z\ = Moreover, the chess¬ 

board estimate holds, see [5j Theorem 4.1]. If, e.g., {/ x } x( =a are operators 
supported on the sites of A 

<n/A<n[wj 1/iAi , (2n 

X eA xGA 

where F x = J^zga /x , and / x z> is a copy of / x attached to the site z if z and 
x belong to the same colored sublattice on the chessboard. If they belong to 
differently colored sublattices, then is a copy of 6f x attached to the site 
z. Eq. (12.71) is obtained by repeatedly reflecting the left side about planes not 
passing through the sites. 

Similar considerations apply to the case in which we have to evaluate 
the average of a product of operators supported on the occupied bonds of 
one of the four classical (columnar) ground states. In this case we repeatedly 
reflect about planes that pass neither through sites, nor through the occupied 
bonds. For these kind of reflections, we need L to be divisible by 4. In the 
next section we will apply this analogue of (12.7p to estimate the probability 


5 











of any Peierls contour in terms of that of the ‘universal Peierls contour’ that 
is as large as A itself, and is much easier to evaluate. 

Note that (12.61) and (12.71) remain valid as A —* oo, in which case the 
average (-) A reduces to the Gibbs measure of the original quantum dimer 
model, simply denoted by (•). From now on we restrict ourselves to this 
limiting case, and we implicitly assume that no unphysical configuration of 
half-dimers appears in the states under consideration. We also drop the oo 
label from H™ and identify H\ with —U — eT + zM, as in (12.lj) (recall that 
we fixed /i = 1). 

2.2 The Peierls argument 

We introduce on-site projectors associated with the 4 classical (e = 0) ground 
states, labelled 1,2, 3,4, as follows: 

if xi is even 
if x\ is odd 

lx if #2 is even 
lx if x 2 is odd 

We want to prove Long Range Order (LRO), in the sense that 

4 

E dib) > !/4, (2.10) 

i=1 

for all A, x, y. The condition (I2.10|) guarantees that the infinite volume Gibbs 
state (■) Z 2 = lim|A|_Kx> (•) is not pure, i.e., it does not satisfy the clustering 
property!: therefore, (12. 10)1 implies that the model exhibits a phase transi¬ 
tion, in that it admits multiple pure states (at least four, by the translational 
and rotational symmetries). To see that (I2.10p violates the clustering prop¬ 
erty, note that (l x ) = (1— Po)/4, Vx e A, Vi = 1, 2, 3,4, where p 0 = (l*). The 
clustering property would imply limix-y^oo JT =1 (ixly) = (1 — Po) 1 2 /4 < 
1/4, which is in contradiction with (12.101) . 

1 The clustering property is the condition that limix-yi^jx, (A x B y ) z2 = (A x ) z2 (B y ) z2 

for any two observables A x and B y supported around x and y, respectively. 



1 * = 


lx = 


lx if X\ is even 
lx if x\ is odd 

if if x 2 is even 
lx if x 2 is odd 


( 2 . 8 ) 

(2.9) 
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Recalling that J2tj=o (ix^y) = 1, a sufficient condition for (I2.10|) is that 
for small enough temperature and monomer density, 

( 1 i 1 y) < S ' I = S' if 0 <i,i <4 and i (2.11) 

< <1») < i Vi = 0,... ,4, (2.12) 

for all A, x, y. The number 21 = 25 — 4 is the number of possible ‘bad pairs’ 
of indices attached to the projectors in x and y. 

In order to bound the left side of (12.111) . we insert at each site z dif¬ 
ferent from x and y a resolution of the identity / = 1 *, and then we 

expand the product of sums, thus getting a sum over configurations of pro¬ 
jectors (we repeat that, since we let A —>■ oo, only physical configurations are 
retained). These configurations are the quantum analogues of the classical 
dimer configurations in m In analogy with the classical case, we associate 
each configuration with a set of Peierls contours mmm, which are closed 
paths on the dual lattice, separating either good regions of different types 
(a good region being a connected region covered by projectors all with the 
same index i £ {1, 2, 3,4}) or a good region from a bad region. We identify 
the contour enclosing the good region containing x and call it 7 . Every bond 
crossing 7 comes with two projectors associated with a bad pair of indices. 
These bonds can be of 4 different types, depending on whether they are of 
the form (x, x + eq) with aq even, (x, x + ej.) with x\ odd, (x, x + e 2 ) with x 2 
even, or (x, x + e 2 ) with x 2 odd. We single out the most numerous among 
these four types, which consist, therefore, of at least | 7|/4 elements; let T 7 
be the corresponding set of bonds. 

For the purpose of an upper bound, we retain only the projectors on the 
bonds in T 7 and throw away (i.e., bound by 1 ) all the others. As a result, 

£ ( n (2-13) 

7>- x {(iz iV)}(z,z')6T 7 (z,z')£T 7 

where the first sum runs over contours enclosing x and the second over the 
bad pairs of indices associated with the bonds (z, z') in T 7 ; if z is the site 
belonging to the same good region as x, then i z = i. 

I 11 order to evaluate (II(zzper ^) > we use ^ ie chessboard estimate 
P, Theorem 4.1], in the way explained at the end of previous section, thus 
finding 

( n p t p y)<d Ti \ ( 2 . 14 ) 

(z,z')eT 7 
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where 


(2.15) 


and 


<5 |a|/2 = max 


2 = max{(P 34 ), (P 32 ), (P 30 ), (P 20 )] 

p‘ s = II U',...i 



(2.16) 


x: rriG4Z 


is a universal projector , in the sense of [5], onto a particular periodic classical 
state. The cases of interest, that is P 34 , P 32 , P 30 and P 20 , are described in 
Figure 1. Each one involves all the vertices in A, not just those in the bonds 
of P 7 . 



Figure 1: A pictorial representation of the four universal projectors that 
are obtained after repeated reflections of the bad bonds crossing the con¬ 
tours. See (12.141) et seq. In the third and fourth pictures the dots represent 


monomers. 












Similarly, the left side of (12. 12p can be bounded by (P 0 ) 1 ^', with P° = 
IIxgA ^x- For later use, we also define P l = Hxga ^x- 
We will prove that 

<P°) 1/|A| < e -ftz+ 1 / 2 ) ? (2.1.7) 

which implies (\2.12\i for z > —1/2 and /3 large enough. Moreover, 

5 < maxje-W 8 - 6 ) + ) 1/12 , e - 0 (*+ 4 / 2 )} (2.18) 

which is as small as desired, provided z > — 1 / 2 , and £ and /3 _1 are sufficiently 
small. 

Assuming this, the validity of (12. lip and, therefore, existence of LRO, 
follows from the standard Peierls estimate: by (j2.13j) . (12.141) . we get that if 
i,j > 0 and i ^ j, and 5 is small enough, 

< v (45) | 7 I /4 (2.19) 

7 >- x 

Now, the number of contours of length 2 n (note that the length of each 
contour is even) enclosing x is smaller than n 2 9 n /12, which gives < 

~bYh n > 2 ™ 2 9 n (4<5) n / 2 . Requiring this to be < 1/28 gives (j2.2j) . provided that 
(12.171) and (j2. 18j) are valid. Our next goal is to prove these estimates! 


2.3 The universal projector 

We first illustrate the estimate for (P 32 ), the ones for (P 34 ), ..., (P°) being 
very similar, or simpler. We proceed as in p, Section I.E and III]. Let 
{V^}i>o be an orthonormal set of eigenfunctions of H\, numbered in order 
of increasing energy (i.e., Eq is the ground state energy). Moreover, let 
(j) be the classical state corresponding to the universal projection P 32 , i.e., 
P 32 = \4>){(f>\- We write 



1 

~Zk 




1 


e fiEi \(il>i\<i>)\ 2 + 

Ei<iE/Q-\-A\A\ i: 


( 2 . 20 ) 

e^ Ei \&i\<t>)\ 2 

®i>®o+A|A| 


R— + P+, 
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where A is a cutoff, which can be fixed to be A = 1/16: the criterion for 
the choice of A is that the energy E 0 + A|A| should be approximately in the 
middle between the ground state energy E 0 and the energy of </>. Neglecting 
£, the ground state energy would be — |A|/ 2 , while the energy of </> would be 
—3|A|/8, i.e., it would be separated from the ground state energy by a gap 
2 A|A| with A = 1/16. 

2.3.1 Estimate of R + 

We start with the easier term, R + . Using the condition E % > E 0 + A | A | we 
get 

R + < J_ e -/3( E o+A|A|)_ ( 2 . 21 ) 

Z\ 

To get a bound on Z\, we restrict the trace to one of the (monomer-free) 
classical ground states, D 0 . Thus, 

Za > (Do I e /3(C7+eT_zM) | D 0 ) > e P(Do\(U+eT-zM)\Do) = & pUo (2.22) 

with Uq = |A|/2. Moreover, T < U, so that Ha > —U(l + e) + zM and, 
therefore, E 0 > — Uq(1 + e). All in all we find 

R + < e _ / 3 ( A l A l _e( 7 o) — e —/3I A I( 1 / 1 6—e/2) ^ (2.23) 


2.3.2 Estimate of 

We now turn to the more subtle R_ term. Here we use the exponential 
localization strategy of [ 6 ]. We bound 

A-< max |(^|</>) | 2 . (2.24) 

The goal is to show that the overlap between the two states is exponentially 
small. It is clear that if e were zero then cj) and would be orthogonal 
because they would belong to subspaces with very different energies. If e ^ 0 
is small, as is the case here, the two states are eigenfunctions of different but 
close Hamiltonians (with and without eT) with very different energies. The 
localization principle says that they are still almost orthogonal. 

Let A = —U + Uq + zM and B = —eT. We also need the projector n 
onto the subspace of eigenvectors of A with eigenvalues larger than |A|A|. 
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Note in particular that 0 is i n the range of II. The eigenvalue equation for 
0 , reads 

(A + -B)0i = (Ei + = A 0,, (2.25) 

so that 0, = —(A — A + i A) _ 1 (i? — iA)0,, for some S > 0. Therefore, (0*|0) = 
— (0,|(£> — iA)(A — A + i5) _ 1 0) and, passing to the limit S —* 0, 

(0#) = -(0 i |S(Al-A)- 1 |0). (2.26) 

Not only 0 is in the range of II, but this is also the case for B(A — A) _ 1 0- 
Therefore, in (I2.26P we can freely add some projectors: 

(0*10) = H^|nsn(A - A) - 1 |0). (2.27) 

We can now iterate this, as long as ( B(A — A) _1 ) n 0 is in the range of II, 
which is the case if, e.g., n < A|A|/3. We thus get 

(0,|0) = (-l) n (0i|(ni?n(A - A) - 1 ) n |0). (2.28) 

Now, ||ILBII|| < e|A|, and ||II(A — A) -1 1| < (3A|A|/2 — A) -1 . Moreover, using 
the fact that — T < U , we have H\ = —U + zM — eT < —U( 1 — e) + zM. 
Therefore, if z + 1/2 > e/2, we get E 0 < —U 0 (l — e), so that, 


A — Ei + Uq < E 0 + Uq + A|A| < sUq + A|A| — |A|(e :/2 + A). 


Consequently, ||II(A — A) *|| < 2|A| X (A — e) l . Using these estimates and 
(I2.26P we obtain, choosing n = A|A|/3, 


wm < (00 


A|A|/3 


(2.29) 


which implies the desired exponential bound on R _. 

By inserting the bounds we obtained on R + , R_ into (I2.20p . we get 


(T > 32)^/|A| < g—/3(2A—e) + ( 2g ) 4A / 3 


(2.30) 


with A = 1/16. By proceeding analogously, we find that ( P 34 ) admits exactly 
the same bound. 

The estimate for (P°) is much simpler: in fact the state 0o corresponding 
to the projector P° is an eigenstate of H\ with energy —Uq + [z + 1/2) |A|, 
which together with (I2.22p immediately implies (12. 1 7p . 
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Finally, we have to estimate (P 30 ) and (P 20 ). By using RP again we can 
bound both of them by (P 0 ) 1 ^. This is accomplished by repeated reflections 
about the vertical planes separating regions occupied by dimers from regions 
occupied by monomers, see Fig.l. In this way, we obtain the inequalities 
(P 30 ) < (P 3 ) 1/2 (P°) 1/2 and (P 20 ) < (P 2 ) 1/2 (P°) 1/2 , but (P 3 ) and (P 2 ) are 
smaller than 1, obviously. Putting things together we obtain the desired 
estimate (12.181) . which concludes the proof of our main result. 


3 Possible extensions and conclusion 

Our model, consisting of dimers with (attractive) parallel dimer-dimer inter¬ 
action and with a ‘flipping’ type ‘kinetic energy’, has been shown to have 
LRO at low temperature and not too large flipping rate. We can also allow 
monomers, as long as their chemical potential is not too negative. 

Unlike other proofs of these properties, ours uses reflection positivity 
(RP), which has the advantage of simplification and the possibility of easily 
obtaining estimates of the allowed constants. It also has the advantage that 
the existence of RP allows one to make statements about the signs of some 
correlation functions and inequalities among them. We do not go into details 
about these matters in this paper, however. In the interest of simplicity we 
concentrated only on the existence of LRO. 

Our work can easily be extended in several ways. One is to other dimen¬ 
sions. The dimers can be replaced by plaquettes in 3 dimensions, although 
they can also remain as dimers. 

Another extension is to give a certain amount of dynamics to the monomers 
by adding a term to the Hamiltonian that annihilates a pair of adjacent 
monomers and replaces them with a dimer (and the reverse, of course). Other 
kinetic terms allowing, e.g., sliding of dimers over unoccupied sites, are not 
generally permitted because they break reflection positivity. 

Extensions to other lattices are possible, too: for instance, the quantum 
dimer model on the hexagonal lattice is reflection positive. However, the 
proof of the chessboard estimate in that case is not as immediate as on the 
square lattice; see [4] for a use of repeated reflections on the hexagonal lattice 
in the context of the Peierls instability in graphene. 
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